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Abstract 

We argue that one can search for physics beyond the standard model through mea- 
surements of the isospin-violating quantity A~° = T(B~ — > p~j) /2T(B° — > p°j) — 1, 
its charge conjugate A +0 , and direct CP violation in the partial decay rates of — > 
p ±r y. We illustrate this by working out theoretical profiles of the charge-conjugate 
averaged ratio A = |(A +0 + A~°) and the CP asymmetry Acp(B — > p ± 7) in 
the standard model and in some variants of the minimal supersymmetric standard 
model. We find that chargino contributions in the large tan (3 region may modify 
the magnitudes and flip the signs of A and Acp(B ± — > p ± 7) compared to their 
standard-model values, providing an unmistakeable signature of supersymmetry. 



Measurements of the radiative decays B — > K*^ H and B X s ^ |2| have trig- 
gered a large number of theoretical studies whose aim is to provide precision tests of 
the flavor sector in the standard model (SM), and to search for possible hints of new 
physics, particularly supersymmetry ||. The related Cabibbo-suppressed decays 
B — > fry, B — > wy and B — > X^J, for which experiments have so far provided only 
upper bounds j|, but which surely will be measured at i?-factories, have also been 
studied at great length. Within the SM, these latter decays are particularly interest- 
ing because they potentially allow us to determine the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix element Vtd, or, more generally, the quark mixing parameters p and 
fj of the Wolfenstein parametrization of the CKM matrix ||. While the inclusive 
decay is theoretically more robust 0, it is experimentally very challenging. In view 
of this, considerable effort has gone into consolidating the theoretical profile of the 
exclusive decays B -> Vj (V = K*, p, w) in the SM [0, g g 0, 0, 0, 

In this letter, we argue that the interference of the short-distance (SD) penguin 
amplitude and long-distance (LD) tree amplitude in exclusive radiative S-decays, 
which is often considered as an impediment to a precise determination of the CKM 
parameters from their branching ratios, may turn out to be a boon in disguise in 
searching for new physics. To illustrate this point, we focus on the decays B°(B°) — > 
p°7 and B ± — > p ± 7, and consider the isospin-violating ratio defined as 

= T(B- -. p- 7 ) _ 

- 2 r(B»^ P » 1 ) ' w 

along with its charge conjugate A +0 . (Since theoretical estimates give t{B ± ) = 
t(B°), to within a couple of percent, and the present data support this conclu- 
sion ||14|| , the quantities A ±0 can be interpreted in terms of the branching ratios.) 
Note that the ratios A ±0 deviate from zero (their isospin limit) due to the SD-LD 
interference effects mentioned above. 

We compare the profiles of the charge-conjugate averaged ratio A = |(A +0 + 
A -0 ) in the SM and in a class of variants of the minimal supersymmetric standard 
model (MSSM) in which all non-diagonal flavor transitions take place essentially 
via the CKM quark mixing matrix. Although the SM and the MSSM yield similar 
values of A in some regions of parameter space, in other regions the MSSM may 
change this ratio significantly. This can happen in two different ways. First, in some 
MSSM models a larger value of the angle a in the unitarity triangle is preferred |15 



Since A increases with a in the quadrant 7r/2 < a < tt P, |TD[, the ratio A may be 
enhanced in the MSSM. The second effect, which is particularly striking, is that the 
sign of A can be flipped in MSSM models. This can happen in that region of large 
tan f3 supersymmetric parameter space in which the chargino-stop contributions are 
known to flip the sign of the effective matrix elements of the electromagnetic and 
chromomagnetic penguin operators [IB, |T7|, |I8, . 



Finally, we also consider the direct CP asymmetry Ac? in the decays B^ p ± 7- 
(The direct CP asymmetries in B° (W) -> p°7 and B°(W) — ► co>7 are very similar to 
^4cp(-B ± — y P ± 7) ) though their time evolution will be modulated by B°-B° mixing 
effects.) We find that, for MSSM's with large tan/3, the sign of Aqp may turn out 
to be opposite that of the SM. 



1 



Having summarized our main results above, we now present the calculation. To 
compute the radiative weak transitions (b —>■ dj), we use the effective Hamiltonian 



n 



(2) 



Here, = V g bV* ql are the CKM factors, and we have restricted ourselves to those 
contributions which will be important in what follows. The operators 0\(p) and 
02 (aO are the four-quark operators 
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2 = {d a T^u a ){u p T ^bp) 
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where Y ^ = 7 M (1 — 75), a and (3 are the SU(3) color indices, and C\ and C 2 are the 
corresponding Wilson coefficients. 07 is the magnetic moment operator 



07 



^(1-75)^6, 



(4) 



where is the electromagnetic field strength tensor. We note that the coefficient 
C?^ (p) also includes the effect of the four-quark operators 05 and 06, and the 
operator matrix elements and their coefficients are calculated at the 6-quark mass 
scale p = m b . For further details and definitions, see Ref. §|. 
The decay amplitudes of interest can be written in the form: 



M(B~ -f p-j) 
M(B° — > p°j) 
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where isospin symmetry has been used in writing a P 



(-) 



,(°) 



= ap for the penguin 



amplitudes, and a is one of the angles of the CKM unitarity triangle. The dynamical 
quantities R L 

are the ratios of the reduced LD and SD amplitudes in the decays B~ — > p~7 and 



and Rl\ which are in general complex due to strong interactions, 



A, 



respectively. In general, ap, ^ and R^' are all model-dependent. 



(0) 



Light-cone QCD sum rules, which take into account the dominant W -annihilation 



and H /± -exchange contributions, typically yield R^ L ^ ~ —0.3 ± 0.07 and R 
0.03 ± 0.01 g g. These estimates, which are obtained using the factorization 
approximation, have been essentially confirmed by a recent calculation in which 
non-factorizable corrections are proven to vanish in the chiral limit to leading twist, 

13]. In addition, long-distance contributions from other 



>(0) 
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in the heavy quark limit 

topologies have been estimated systematically and found to be small 
Eventually, radiative decays B ± i ± v^ can be used to compute the leading {W ± - 
exchange) topologies in a model- independent way |0| . Of course, one still needs to 
know ap to get the branching ratios. 

The expression for the ratio of the branching ratios of interest can be written as 



B(B° p° 7 ) 



1 H tti r" 



eff 



A 



(rf) 
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where r p u lumps together the dominant (W-annihilation) and possible sub-dominant 
LD contributions. Borrowing the notation from Ref. fl3| . 



ATi 2 m p ±C 2 + C l /N c 



eAe i4> A = - ',;;r > r 



m b 



eff 
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and noting that A,[fYA 



{d) 



*, which holds in the SM and in the 
MSSM models being considered here, the isospin breaking ratios [Eq. ([!])] can be 
expressed as 



A 



±o 



2e t 



cos (f> A F 1 =f sin0 A F 2 + ~e A {F? + F 2 2 ) 



Here, F 12 are (implicit) functions of the Wolfenstein parameters p and r\: 



AW 



cos a 



F, = - 



sin a 



(9) 



Ad) 



with (F 2 + F^ = AW/A| 

The charge-conjugated averaged ratio, defined as A 
following leading-order (LO) expression: 
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A +0 1, has the 
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where the near equality reflects that, in this approximation, the strong interaction 



phase <f> A disappears in the chiral limit [13 



In fact, one can go to next-to-leading-order (NLO) in the calculation of the above 
quantities. The NLO-corrected expression for the branching ratios and A can be 
derived from the corresponding calculations for the inclusive decay B — ► X s 7 [2fJ |21 
and B^X dl M: 
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Here Gp is the Fermi coupling constant, a = a(0) = 1/137, Tf is the B 
form factor involving the magnetic moment operator Oy, evaluated at q 2 = 0, and 
■(°) e //_ ^he quantities A^j and 7 represent the real and imaginary parts 



of the explicit 0(a s ) contributions to the matrix elements evaluated at a scale \i\ 
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where r^s are complex numbers. Expressions for the various quantities appearing 
in the above equations can be found in Refs. |20|, 2l| . We stress that the gluon 



bremsstrahlung parts have been dropped in calculating T(B — > fry), except those 
needed to cancel the divergence in the 0(a s ) virtual corrections in the decay b — > dj. 
Note that, in the above rate, all terms higher than 0(a s ) have to be dropped for 
theoretical consistency. The expression for T(B° — > p°7) can be obtained by obvious 
replacements, except that L R (B°) <C L^(_B ± ). 

Using the above expression, the NLO isospin-violating ratio A is found to be: 



Anlo — Alo 
2e A 



(0)e// 



FiA<S* - (F| - Ff) Al + e A (Ff + F^) ( A ™ + F^)] (14) 



where Alo is given in Eq. (|i~0|). 

The values for the various input quantities used in the numerical calculations 
of Alo and A NLO in the SM are as follows: df )eff (m b ) = -0.318, A^ ]t = -0.022, 
A 1 ^ = +0.049, and e A = —0.3. The remaining ingredient is a determination of 
the allowed ranges for the functions Fi and F 2 . Taking into account the present 
experimental and theoretical constraints on the parameters of the CKM matrix, the 



profile of the unitarity triangle in the SM was presented by two of us in Ref. |Tj|. In 
Fig. [I] we show the allowed iq-a and F 2 -a correlations at 95% C.L. In these figures, 
the SM plots are found on the left-hand side, and are labelled by / = 0.0. The 



ranges of the hadronic parameters fB d \J B Bd and B K used in these fits are indicated 



on top of the figures. (For definitions and further discussions, see Ref. Note 



that the CP phase a is constrained to lie in the range 75° < a < 121° at 95% C.L. 



15 



With this information, we can now calculate the ratios Alo an d Anlo m the SM. 
In Fig. ^| the results are shown for these quantities as a function of the angle a. In 
these figures we have assumed that |V^b/Vtd| = 0.48 (its central value ]I5| ). However, 
for a given value of a, Alo an d Anlo may in fact take a range of values. This 
residual CKM-related range is given essentially by the iq-a correlation presented 
in the upper-left plot in Fig. [1|. Note that the isospin-violating ratio is very stable 
against NLO corrections in the SM. This observation, together with the discussions 
earlier about the determination of e A , makes A suitable for precision tests of the 
SM. In particular, its measurement will determine a in the SM. 

We now turn to the direct CP asymmetry WM. As noted earlier, the strong 



interaction phase (f>A of Eq. (^) disappears in the chiral limit |13| , which implies 
that, to lowest order, there is no CP-violation in the decay rates for B — > /ry. 
Therefore, the strong phases in the exclusive decays B ± — > p ±r y and B°(B°) — > p°7, 
which are necessary for inducing direct CP-violation, must be generated by higher- 
order perturbative QCD corrections. Concentrating on the charged B decays, we 
define the CP asymmetry as 

_ B{B~ -> pjy) - B{B+ -> p+ 7 ) 
CP ~ B(B~ -> p~i) + B(B+ -> p+ 7 ) ' 1 ' 
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Since, in the heavy quark limit, there are no non-factorizing strong phases in the 
jy-annihilation part of the B — > amplitudes ||13|| , the strong phases are generated 
entirely by the Bander-Silverman-Soni mechanism p3[, which involves the interfer- 
ence of the penguin operator O7 
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This 



and the four-quark operator C 2 
mechanism has been employed by Greub, Simma and Wyler to calculate *4.cp, using 
a wave function model |25j for the mesons. Since we are working to leading twist, we 
shall ignore the effects involving virtual corrections off the spectator quarks, arguing 
that they are suppressed by powers of l/m&. In that case, the CP asymmetry is 
determined by perturbation theory up to a non-perturbative quantity which can be 
determined from the ratio A. The expression for Aqp is given by 
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2Fo 



CP 
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The quantities A^* and Aj take the values AY'" = -0.016 in the SM and Af = 
+0.046 in both the SM and MSSM. 

Note that A and ^4cp are complementary measurements. Dropping the small 
0(e\) terms in Eq. (|14|), we see that A is essentially proportional to Fi, while Aqp is 
proportional to F 2 . Thus, for a ~ tt/2, A is very small, while Acp takes its maximal 
value. Conversely, if the value of a is far from tt/2, the CP asymmetry decreases, 
while A becomes measurable. 
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Figure 1: Upper-left: SM F\-a correlation. Upper-right: MSSM Fi~a correlation. 
Lower-left: SM F 2 -a correlation. Lower-right: MSSM F 2 -a correlation. The pa- 
rameters used in calculating the correlations are indicated on top of the figures. 
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Figure 2: The isospin-violating ratio A at LO and NLO in the SM. We have set 
€a = —0.3. The curves correspond to the central value of the CKM fits IV^/V^I = 
0.48 pi. 



We are now ready to examine the supersymmetric contributions to A and *4cp- 
To begin with, we note that the NLO corrections to the decays B — > Xgj have 
been calculated in only one particular realization of the MSSM, the so-called min- 
imal flavor violation scenario [[T7|]. While this calculation considers an important 
parameter space in the MSSM, it nevertheless neglects other contributions, such as 
those from gluinos, which are important in other regions of parameter space [I8|. 



In the small-tan /3 domain, where the neglected contributions are small, we have 
numerically calculated the NLO quantities and found that the NLO correction to 
A in the MSSM with minimal flavor violation is very similar to that in the SM, 
and hence unimportant. The complete NLO corrections for the large-tan (3 case, 
including gluino contributions, are not yet available. Hence, in comparing the SM 
profile with that of the MSSM, we shall restrict ourselves to Alo- 

Supersymmetry can affect A and ^4cp i n two distinct ways. First, the allowed 
values of the functions Fi and F 2 are different in the MSSM. We recall that the 
supersymmetric contributions to the mass differences M 12 (-B) and M\ 2 (K) can be 
written as follows (for details and references, see Ref. [|T5] ): 



AM d = AM d {SM)[l+f d {m x ±,m iR ,m H ±MnP)), 
AM, = AM s (SM)[l + f s (m x ±,m iR ,m H ±,t&nf3)}, 

+ rjttyth(yt)[l + f e (m x ±,m~ t2 ,m H ±,tanf3)}A 2 \ 4 (l - pfj . (17) 

To an excellent approximation, one has fd = fs = fe = f- The quantity / is 
a function of the masses of the (lighter) right-handed top squark (mi ), chargino 
(m-±) and the charged Higgs (m#±), as well as of tan (3. The maximum allowed value 
of / depends on the model. Typical values are: minimal supergravity (/ = 0.2), 
non-minimal supergravity (/ = 0.4) (|TJJ], and MSSM with constraints from electric 
dipole moments (EDM's) (/ = 0.6) f26|. The plots in the upper right-hand and 
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lower right-hand corners of Fig. [I] show the allowed Fx— a and F^-ol correlations, 
respectively, for the MSSM with / = 0.6. We see that these correlations can be 
measurably different from the SM. In particular, much larger values of a are allowed 
compared to the SM. Thus, for / = 0.6, the fits yield 86° < a < 141° at 95% C.L. 

The second way in which supersymmetry affects A and Ac? is via the Wilson 
coefficients. In contrasting the SM and MSSM profiles, we assume, as per the 
usual expectations, that the coefficients of the tree amplitudes, Ci and are 
the same in these models, but that C\ ? (p) may differ. This latter coefficient 
is constrained by the measured branching ratio of the decay B — > X s 7, yielding 
a bound 2.0 x 10" 4 < B(B -> X sl ) < 4.5 x 10" 4 at 95% C.L. §. The resulting 
constraints on the magnitude and phase of the ratio C^ e ^ jQ^ff^w) j n ^ e context 
of the MSSM being considered can be summarized as follows. In the absence of the 
constraints on the EDM's of the neutron and electron, the real and imaginary parts 
of Cj^/C^ S can vary substantially. But if one takes into account the EDM 
constraints, the imaginary part of cf^ e ^ j(^ e ff^ M > j s highly suppressed. However, 
depending on the value of tan /3, both positive- and negative-valued solutions of C 7 
are allowed. For example, a recent analysis of C^ e ^ /cj ) e ^( SM ) j n the minimal 
supergravity model yields values in the ranges 0.7 < RefC^ 6 '^ /C^ S ] < 1.2 for 
small tan/3 (say tan/3 < 10), but for larger values of tan/3, negative values of this 
ratio are admissible. Thus, for tan/3 = 30, a range —1.5 < Cf^ e ^ ^ / Cf^^ '^ SM ^ < —0.8 
is allowed by present data |19| . 



As for Acp(B —> p 7), the dominant contribution, proportional to AJ, is iden- 
tical for the SM and MSSM. However, the value of the quantity A^p 1 in the MSSM 
depends on the region of parameter space considered. For small tan/3, A^p 1 has 
almost the same value as in the SM, while for large tan /?, it may appreciably differ 
from the SM value. However, since its contribution to AcviB^ — > p ±r y) is sup- 
pressed due to the factor, its precise value is not so important numerically. In 
calculating Acp in the MSSM, we have set A^ 1 = 0. 

In Fig. |3] we contrast the expectations for A in the SM and in two variants of the 
MSSM, characterized by small and large values of tan /3. The residual CKM-related 
range in the allowed values of A is again given essentially by the F\—ct correlation 
presented earlier in the upper two plots of Fig. [I] for the SM (/ = 0.0) and MSSM 
(/ = 0.6). Note that if the large tan/3 MSSM solution is realized in nature, the 
measured value of A ±0 can be markedly different than in the SM. More importantly, 
since jC^ 1 S can be negative, the sign of A ±0 may change over a large 

region of the allowed CKM parameter space. This would be a striking signature of 
new physics, and would strongly suggest the presence of supersymmetry. 

Similar effects appear in the CP asymmetry Ac?. Since, in the models being 
considered here, there are no other phases at this order, Eq. holds in the 
MSSM, with the proviso that numerically C^ e ^ now depends on the parameters 
of the MSSM. In particular, for large tan /3, the CP asymmetry can be significantly 
larger than the one in the SM. Again, since Cf )eff (MSSM) ~ -Cf )eff (SM) is 
allowed, the CP asymmetry Acv(B ± — > p ±r y) reverses sign in this case (as does the 
asymmetry in the inclusive decays ^4cp(-B ± - ► 7)). This is illustrated in Fig. BL 
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Figure 3: The isospin- violating ratio A in LO (left) and Acp^B^ — > p ±r y) 
(right) with e A = -0.3 in the SM (solid line), and in the MSSM with 
(tan/5,Cf e// /Cf e//(SM) ) = (3,0.95) (dot-dashed line), (tan/5, 4° )e// /4° )e//(SM) ) = 
(30,0.8) (dashed line) and (tan/?, C f )eff /cf )effm) ) = (30,-1.2) (dotted line). The 
SM and MSSM curves correspond respectively to \V u b/V t d\ = 0.48 and \V u b/V t d\ = 
0.63, which are the central values in the CKM fits [|TjJ. The allowed ranges for a in 
the MSSM from these fits for / = 0.6 are also indicated. 



Once again, this would be a clear signal of supersymmetry in the large tan (3 domain. 

In summary, we have examined the effects of supersymmetry on two observables 
of B — » p7 decays: the isospin-violating quantity A and the direct CP asymmetry 
^4 CP . We find that, in MSSM models, the predictions for these quantities can be 
substantially modified. In particular, over a large region of parameter space, the 
signs of A and Acp may be flipped, which would be a clear signal of new physics, 
and would point directly to the presence of supersymmetry. 

We thank Christoph Greub, Toru Goto, Gudrun Hiller, Klaus Honscheid, David 
Jaffe, Alexey Petrov and Dan Pirjol for helpful discussions and communication. The 
work of D.L. was financially supported by NSERC of Canada. 
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